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The arrow of time and the accelerated expansion are two fundamental empirical facts of the Uni¬ 
verse. We advance the viewpoint that the dark energy (positive cosmological constant) accelerating 
the expansion of the Universe also supports the time asymmetry. It is related to the decay of meta¬ 
stable states under generic perturbations, as we show on example of a microcanonical ensemble. 
These states will not be meta-stable without dark energy. The latter also ensures a hyperbolic mo¬ 
tion leading to dynamic entropy production with the rate determined by the cosmological constant. 


I. INTRODUCTION 

Astronomical data point out the existence of a positive 
cosmological constant A (dark energy). 

The precise origin of A is not yet clear. Various mod¬ 
els attribute it to macroscopic vacuum fluid, geometrical 
term, scalar fields or modified gravity [l|. An influential 
scenario for A > 0 is that it emerges due to vacuum fluc¬ 
tuations which are able to induce negative pressure Q; 
see e.g. ii for the numerical fit of estimations and the 
observed density of dark energy. The positive A is among 
the necessary conditions along with the 2nd law also in 
the Conformal cyclic cosmology d, @1 ■ 

The time is ripe for asking about the implications of 
A > 0 in basic physics. We aim to show that A > 0, 
in contrast to A < 0, leads to a specific mechanism for 
an emergent thermodynamic arrow of time and entropy 
generation. It is related to parametric instability of the 
bound, gravitating motion that becomes metastable with 
respect to generic perturbations in the presence of dark 
energy. 

The time asymmetry and the 2nd law of thermody¬ 
namics are among the basic empirical facts on the Na¬ 
ture. It is known that laws of physics are invariant with 
respect to time-inversion [7Hl5l| (more precisely they are 
only CPT-invariant; the difference between T and CPT 
is an important one [l^ . but it is relevant only for high 
energies [13). However, there are several fundamental 
types of motion whose time-inversion is not observed. 
They are called arrows of time [M3: 

- thermodynamical: increase of entropy in a closed 
evolving system; 

- electrodynamical: physics is dominated by retarded 
potentials, although advanced potentials are formally al¬ 
lowed, they are not observed; 

- cosmological: expansion of the Universe; 

- quantum-mechanical: the apearance of definite mea¬ 
surement results accompanied by reduction of quantum 
state. If quantum mechanics is an emergent theory, this 
arrow may be caused by a sub-quantum one; see [22, [2^ 
for possible scenarios. 

There are certain relations between the arrows [ 13 ; 
in particular, the quantum-mechanical arrow can to an 


extent be reduced to the thermodynamic one [T3, fl3 - 
Recent research clarified the place of this arrow in mi¬ 
croscopic dynamics [Sill and its relation with external 
perturbations [2l|. 

In all arrows there are two closely related aspects: ini¬ 
tial conditions and the proper dynamics. Let us recall 
and illustrate this point via the emergence of the ther¬ 
modynamic arrow within the system-bath approach m 
from the T-invariant Hamiltonian dynamics; se e Illl - ll4 
for a general background. It was argued in that 
(i) the thermodynamical time arrow in the system can 
arise in the system due to the limited observability of 
the bath; (ii) while the initial conditions are necessary for 
the emerging of a pre-arrow, the full time arrow will be 
established if also the dynamics of the system is Marko¬ 
vian (no-memory). Namely, when a quantum system S 
interacts with a thermal bath B, the total Hamiltonian 
is split H = Hs -\- Hb + Hi, between, respectively, S, B 
and interaction. The state of S -|- B is described by the 
density matrix D(<) and the von Neumann equation 

ihdtV{t) = HV{t) - V{t)H. 

When the initial state at t = 0 can be split as 
D(0) = Ds(0)Gi?B(0), 

and the state of the system at arbitrary time t is given by 
partial density matrix Ds{t) = trBD(t), where trB is the 
trace over the Hilbert space of the bath, one can see the 
emergence of the thermodynamical arrow in the system 
due to the bath’s incomplete observability. 

In [23 we also discussed the hyperbolicity as a pos¬ 
sible mechanism for the Markovian dynamics. One sce¬ 
nario for this relates to the voids—underdense regions 
in the Universe—that are able to induce hyperbolicity of 
the null geodesics even if the global spatial curvature is 
zero (i.e. in the flat Universe) (HL The properties of 
the cosmic microwave background 26| appear to fit the 
observed void structure in the large scale galaxy distribu- 
tion, iiicluding in the case of the Gold Spot as a supervoid 

[23I23. 

We now make the next step in that approach of the 
emergence of the time arrow i.e. involving one more fun¬ 
damental empirical fact, the dark energy. 
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We adopt an Ansatzth&t the dark energy acts as a bath 
for the observed Universe, thus supporting the emergence 
of the time arrow. The system-bath interaction has to be 
small in order not to distort the state of the system. This 
agrees with the empirical situation of the dark energy 
when the role of its influence on typical macro-physical 
processes remains unnoticed both due to the low value of 
its energy density and the still ambiguous cross-section 
of interaction with elementary particles. In this sense 
our laboratory physics reveals itself within intermediate 
scales on which both the cosmological constant and vac¬ 
uum fluctuations are not easily noticed, although possi¬ 
bly being themselves mutually linked. 

Here the relation between the arrow of time and dark 
energy (A > 0) will be established for a particular 
scheme, though various considerations of this Anzats can 
be possible. Within this scheme, the dark energy facili¬ 
tates the thermodynamic arrow of time on those scales, 
where its influence can be comparable to gravity. 


II. SET-UP 


We consider the limit of the Newtonian gravity, where 
N non-relativistic test particles move in a field of a large 
mass M. Here A reveals itself as an additional parabolic 
potential imposed on the usual inverse-square-law inter¬ 
action [29l - [^ 1^ . As usual in equilibrium statistical me¬ 
chanics, we assume that initially the Wparticle system 
performs a finite motion and can be described by a micro- 
canonical, equilibrium distribution at some fixed energy 
[s^ 1^ . So no arrow of time is present initially. Our aim 
is to look for two scenarios of perturbing this system such 
that for A < 0 (negative or zero cosmological constant) 
the system is stable and continues to perform a finite mo¬ 
tion. These scenarios amount, respectively, to fast and 
slow perturbations. Moreover, for the slow perturbation 
scenario A < 0 implies that the system returns to exactly 
the same state as before the perturbation. However, for 
both considered perturbation scenarios, A > 0 (positive 
cosmological constant as observed in our universe) leads 
to changing the finite motion to an infinite one so that 
the ergodivity is violated, i.e. the system moves in one 
ergodic component, and the real motion is not anymore 
similar to the time-inverted one. Consequently, the dy¬ 
namic entropy increases with the rate oc \/A. 

Since all these effects relate to A, we shall choose to 
work with the simplest situation N = I: one-particle 
system for which the initial microcanonical distribution 
is well-defined. All the obtained effects exist also for 
N >1. 

Thus for a test mass m in the field of a larger mass 
M :S> m the Newton equations read: 


R = -d^{R)/dR, 


where R is the interparticle distance, and 
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Here m^{R) is the potential energy of the test parti¬ 
cle, L is the (constant) orbital momentum, — is the 
gravitational attraction, and — R^ is the potential 

generated by cosmological term. It is characterized by 

Pv = Ac^/(87rG'), (3) 


the mass density of the vacuum fluid, if A > 0 (dark 
energy) is interpreted in this way. Note that in —^{M -|- 
contribution from the dark energy is seen 
to arise from a homogeneous distribution with density 
Pv- 

The recent estimate for the dark energy density by the 
Planck’s data yields 0.69 fraction of the total density [s^. 

Thus the conserved energy of the test particle is a 


E=^ + $(i?), (4) 

The terms — and — in ([5]) are similar to 

each other [^: they both hold the Newton’s shell theo¬ 
rem (they are the only potentials having this feature) and 
they possess an additional symmetry leading to closed or¬ 
bits. 

Let us mention somewhat different interpretation of 
(011: they apply to a test mass in the homogeneous, 
isotropic universe , where the motion of the test mass 
is influenced only by the matter mass M and the “vac¬ 
uum” mass rA inside of the sphere with the radius 

R. 

Also, the third term in 0 leads to the A-term in Fried¬ 
mann equation [sij: 


■ 2 1/ A, 

a =—{const — na +—), 

CL O 


( 5 ) 


where a is the scale factor and k is a constant. This 
approach to the large scale limit for the Newtonian po¬ 
tential removes the infinities peculiar to the purely New¬ 
tonian cosmology . Its radial dependence can become 
a subject of dedicated astronomical testing based on the 
dynamics of galactic halos, galaxy groups and clusters. 

This is related also to the already discussed scale (e.g. 
0 ), when N-body effects become comparable to the 
dark energy one. The possibility of observin g th e dark 
energy at this scale was recently discussed in [^ . 


III. PERTURBATIONS 

We turn to a detailed investigation of 0. First of all, 
note that the term 2 ^^ in 0 is needed for ensuring a 
bounded motion of the test particle: otherwise, for L = 0 
it will simply fall into the central mass. Apart of that, 
the term does not play any important role in our study 
and for simplicity we replace it by a hard-wall imposed 
at relatively small distance Rq. 

Let us introduce characteristic scales E and R for 
the energy and distance respectively, and write 4>(i?) = 


3 
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E(j){r) in terms of dimensionless r = R/R and 4>{r)'. 

(j){r) =—ajr — (3r‘^ 12, r > rp. (6) 

where a rx M, /S rx pv, and r > tq is the hard-wall 
condition. Fig. [T] shows the form of this potential. It is 
maximal at 


r, = ai/3^-1/3 


= (7) 


Hence the energies ^(rc) > e > ^(ro) [^(?’c) < e] refer to 
bounded [unbounded] motion. 

How the situation changes when the mass M is time- 
dependent? This is the only natural parametric pertur¬ 
bation for this problem. Indeed, if several large masses 
are there (and the test particle moves in the effective field 
generated by them) the inverted harmonic potential act¬ 
ing on the test particle should originate from the inertia 
center of the overall system oc — 32|, and hence it 

cannot become (parametrically) time-dependent (the in¬ 
ertia center is at rest). We stress that the time-dependent 
mass M (t) always stays much larger than the mass of the 
test particle. 

Now provided that the test particle’s motion is 
bounded and M{t) is slow, E{t) can be described via 
the (adiabatic) invariant phase-space volume [H, l40l| 


dRdP'd 


E{t)-^-<^>{R,M{t)) 


( 8 ) 


where P is the momentum [cf. (HJ], and d[a:] = 0 for 
a; < 0 and i?[x] = 1 for a: > 0 is the step-function. The 
conservation of ([8]) relates to the fact that (in the present 
case) the bounded motion is ergodic [4l| . Recall that the 
entropy of a microcanonical equilibrium state is defined 
via the logarithm of ([8]) j^. Its adiabatic conservation 
relates to the second law (4ll| . Eq. m and its generaliza¬ 
tions appear as well in the control theory (d^ . 

Integrating over P in ([8]) and going to the dimension¬ 
less quantities we get from dH]) that [up to irrelevant con¬ 
stants] (|5]) reduces to 


J = 


r(*) ,- 

/ dr \/e- (()(r,a(t)), 

j rn 


(9) 


where e(t) = E{t)/E, and r(t) is the maximal distance 
for the finite motion at time t: 


e{t) = (j){f{t),a{t)). (10) 

Note that f{t) is always smaller than the largest possible 
distance of the bounded motion for a given a{t) and /3(t): 

r{t) < rdt) = (11) 

Thus when changing a from ai to a 2 , the final energy 
£2 is to be determined from 


( 12 ) 


-2 
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FIG. 1: The effective potential versus distance r for ro = 
0.05. Dashed black curve: a = P = 1. Black curve: a = 0.5, 
/3 = 1. The energy e = — 2 refers to a bound motion for a = 
/3 = 1, but when slowly decreasing a{t) this motion becomes 
unbounded, i.e. its energy raises above the red-dashed (upper 
horizontal) line; cf. Fig. [2] 



FIG. 2: The adiabatic invariant J(e) as a function of the 
dimensionless energy e for a = P = 1 (dashed black curve) 
and for a = 0.5, P = 1 (black curve). For both curves pdc) > 
e > (f>{ro), where d>(rc) (</)(ro)) is the largest (smallest) energy 
for the bounded motion. Now energies e £ (-1.5,-2.4889) 
(this interval is denoted by two red, bold points) that refer to 
finite motion under a = p = 1 correspond to infinite motion 
when a slowly decreases from a = ltoQ: = 0.5. The range t G 
(—1.5, —2.4889) is additionally illustrated by the horizontal 
(dashed red) line; cf. Fig. [T] 


A time-dependent a{t) cx M{t) leads to the following 
two scenarios by which the bounded motion can turn to 
unbounded one. Both of them do not exist for A oc pv < 

0 . 

1. Let M{t) oc a{t) decrease slowly. This can model 
slow evaporation taking place from the mass M. It is 
now possible that the bounded motion with sufficiently 
high initial e becomes unbounded; see Fig. [TJ This is 
related to the fact that m does not have solutions £2 
for a range of £1 that initially were sufficiently close to 
4>{rc); see Fig.j^l In other words, during the slow decrease 
of a{t), e{t) grows faster than (j){rc{t)). Now it is crucial 
that the change of a{t) is not very fast; otherwise e{t) will 


J(ai,£i) = J(a2,e2). 
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FIG. 3: The effective potential (j){r) versus distance r for ro = 
0.05 [cf. ((6|]. Dashed black curve: a = /3 = 1. Blue dashed 
line refers to 4>c = —(for a = (1 = 1); cf. 0. For 
all energies below (above) 4>c the motion is bound (unbound). 
Red line: an example of finite motion at energy e = —2. 
Black curve: ((>(r) for a = 2, /3 = 1. When a slowly changes 
from Q = 1 to q: = 2, the energy (given initially by the red 
horizontal line) decreases and always refers to finite motion; 
for Q = 2 it is given by the green (lower horizontal) line and 
is equal to —5.7593, as found from m- But if a changes 
sufficiently fast, the initial energy does not change much 
and now it corresponds to an unbound motion; compare red 
dashed line with the black curve. 
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FIG. 4: Phase-space density for the inverse parabolic poten¬ 
tial with Hamiltonian H(r,p) = where (r,p) are 

canonically conjugate coordinate and momentum; cf. m- 
Green square [(2,4); (2,4)] is the phase-space volume at the 
initial time t — 0. Purple figure is the phase-space volume 
at time t = 3. The initial volume was filled by 10® random 
points; each point was given a small but finite width, defining 
a coarse-graining procedure. The evolution of each point was 
then followed from f = 0 till f = 3. It is seen that due to the 
coarse-graining the phase-space volume increases. 


not change much and will stay bounded. Again pv > 0 
is crucial for this scenario. For A < 0, the influence of a 
cyclic change of a{t) on the system can be made arbitrary 
small, provided that it is sufficiently slow. 

2. Let M{t) oc a{t) increase. This can account for 
situations with accretion driven increase of M. Now 
Tc increases [see ([T])] and e{t) decreases. But if a{t) 


changes suddenly, e(t) will not change much and the mo¬ 
tion will become unbounded; see Fig. [31 (The energy 
e{t) will not change, since during a sudden perturbation 
the force changes by a hnite amount in a small time- 
interval. Hence the acceleration also changes by a finite 
amount, while the changes of the coordinate and veloc¬ 
ity are negligible; this is the general feature of sudden 
perturbations.) 

Thus, there is a larger class of perturbations such 
that a{t) changes not slowly, and e{t) decreases slower 
than (j){rc{t)); see Fig. [31 Hence it is possible that 
e{t) > 4>{rc(t)) at some time t, and the motion becomes 
unbounded. This conclusion may seem counter-intuitive: 
once the mass M increases, the attraction towards the 
center becomes stronger, but the test particle can escape 
the attracting center benefiting from the repulsion by the 
dark energy inside the shell. 

It is crucial for this scenario that a{t) increases not 
slowly. Otherwise, m predicts bounded motion; see 
Fig. [3] This scenario of parametric instability is due to 
Pv oc A > 0, e.g. the situation with A = 0 (no dark 
energy) or with A < 0 (negative cosmological constant) is 
stable with respect to this parametric perturbation. Note 
that when a{t) returns to its initial value—i.e. when the 
perturbation is over—the motion is not turned back to 
bounded. 

Thus in both scenarios the motion changes from 
bounded to unbounded. The latter is not ergodic, e.g., 
only one component of the momentum space is explored. 
(For example, in the Id situation the momentum space 
P has two components P > 0 and P < 0.) Similar ex¬ 
amples of irreversibility generated by non-ergodicity were 
analyzed in (4ll| . 


IV. ENTROPY PRODUCTION 

What happens with the test particle once it escapes 
the bounded part of the potential? The dominant part 
of potential is now —/3r^/2 and the motion generated by 
Hamiltonian H = — f3r^), 

r = p, p = Pr, (13) 

is hyperbolic: it has a positive Lyapunov exponent 
VP oc \/A. Due to the second Lyapunov exponent 
— VP it is phase-space volume preserving as any Hamil¬ 
tonian motion. The exponent VP relates to an expand¬ 
ing eigenspace which is also the stable manifold; see 
Fig. m A small coarse-graining will thus lead to an in¬ 
creasing phase-space volume (and thus growing entropy) 
as demonstrated in Fig. S) It is expected that the rate 
of this increase (i.e. entropy production) will be deter¬ 
mined by VP oc VA. It was shown that under a weak 
noise (which is equivalent to a specific coarse-graining) 
the positive Lyapunov exponent defines the rate of en¬ 
tropy increase for the motion in the inverted parabolic 
potential (d^ . 
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Recall that a coarse-graining (or alternatively an ex¬ 
ternal noise) is standardly needed f or g etting an arrow of 
time for a Hamiltonian dynamics [uHlSj- Hamiltonian 
systems displaying the arrow of time are those, where 
a small coarse-graining leads to entropy increase with 
a rate that only weakly depends on the coarse-graining 

[nHll. 

Note that the oc \/A scaling of the entropy production 
is natural, since it indicates on the inapplicability of the 
whole time-arrow scenario for A < 0. Indeed, for A < 0, 
the cosmological constant amounts to an overall confining 
potential and thus increases the stability of the motion. 

However, this motion is not chaotic (not even ergodic). 
It generically does become chaotic provided that its mo¬ 
tion becomes bounded at some scale; see [3 for a con¬ 
crete scenario related to the inverted parabolic potential. 
It is conceivable that the accelerating particle will meet 
other masses, and the resulting interaction will achieve an 
effectively bounded phase-space. Similar and more com¬ 
plex examples of hyperbolic motion in self-gravitating 
systems were studied in [H, |4^ . There is another argu¬ 
ment showing that the motion in the inverted-parabolic 
potential will change its dynamic regime: for a very large 
R the potential ^{R) will assume very large absolute val¬ 
ues [cf. ([2])], which violates the known applicability con¬ 
dition |$|/c^ <C I of the non-relativistic description (d^ . 

It is known that in a Universe with a positive cosmo¬ 
logical constant A > 0, the long-time evolution of the 
space-time will be dominated by A [d^. Locally (but 
not globally) this Universe will look like empty, since the 
matter will escape through the horizon [^. Now for 
this Universe it was shown that the generalized second 
law holds: the sum of the matter entropy and geomet¬ 


ric entropy (related to the horizon) does not decrease 
(d9l - l5i| . For very late times, where locally (almost) no 
matter is present, the matter entropy will tend to zero, 
while the geometric entropy saturates at a value oc I/A 
(dol - l^ . Our consideration concerns the matter entropy 
and refers to the opposite limit of sufficiently early times, 
where the matter is abundant, there are bound systems 
that perform finite motion etc. 


V. SUMMARY 

Thus, the cosmological constant term induces hyper- 
bolicity of motion with the time asymmetry and dynamic 
entropy production on the large-scale when A is dominat¬ 
ing. We provided a mechanism by which at those scales 
A facilitates the thermodynamic arrow. We stress that 
this mechanism does not directly apply to early universe. 
Indeed, our starting point is a closed (microcanonically 
equilibrium) system that does not show an arrow of time 
before perturbed externally. Such systems exists inde¬ 
pendently on initial conditions of the early universe. 

We also note that the mechanism cannot (and should 
not) explain all occurences of the thermodynamic arrow. 
However, note that even when the dark energy (cosmo¬ 
logical constant) does not dominate the mean density 
(early universe or today’s laboratory scale), it still ex¬ 
ists. To give an example: for a quantum system in a 
laboratory A (= vacuum energy) is not dominating, al¬ 
though it exists (e.g. Casimir effect). Importantly, the 
dark energy can serve as an ideal thermodynamic bath, 
since it does not get any back-reaction. 
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